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ABSTRACT 

The  atomic  structures  of  four  Ca4o+xMg25Cu35.x  (X  =  0,  5,  10,  and  20  at.%)  ternary  metallic 
glasses  have  been  detennined  using  a  synergistic  combination  of  neutron  diffraction,  quantum 
molecular  dynamic  (QMD)  simulation  and  constrained  reverse  Monte  Carlo  modeling.  The 
amorphous  structure  of  these  alloys  can  be  described  as  close-packing  of  efficiently-packed  Cu- 
centered  clusters  that  have  Ca,  Mg  and  Cu  atoms  in  the  first  coordination  shells.  Experimental 
evidence  is  given  for  the  presence  of  Cu  atoms  in  inter-cluster  sites  for  the  two  most  Cu-rich 
glasses.  The  close-packed  arrangement  of  clusters  is  present  only  within  length  scale  a  length 
scale  of  — 10  A,  providing  a  characteristic  medium  range  order  pre-peak  at  Q  ~  1.2  A  in  the  total 
scattering  structure  factors  of  these  alloys.  An  average  coordination  number  of  10  (with  about  5- 
7  Ca,  2-3  Mg  and  1-2  Cu  atoms)  is  most  common  for  the  Cu-centered  clusters.  The  average 
coordination  numbers  around  Mg  and  Ca  are  12-13  (~  6-8  Ca,  3  Mg  and  1-4  Cu)  and  13-15  (7-9 
Ca,  3-4  Mg,  and  2-5  Cu),  respectively,  and  they  are  composition  dependent.  Strong  interaction  of 
Cu  with  Mg  and  Ca,  which  result  in  pair  bond  shortening,  has  been  identified.  Icosahedral  short 
range  order  does  not  dominate  in  these  amorphous  alloys,  although  polytetrahedral  packing  and 
five-fold  bond  configurations  resulting  in  pentagonal  bi-pyramids  have  been  found  to  be  the 
most  common  nearest  atom  configurations. 

PACS  numbers:  6143. Dq,  61.05.fm,  61.25.Mv,  61.43. Bn,  83.10.Rs 
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1.  INTRODUCTION 

It  is  well  established  that  Ca-Mg-  based  metallic  glasses  have  very  good  glass  forming  ability 
(GFA)  [1,2, 3, 4, 5, 6].  Based  on  two  simple  metal  elements,  these  glasses  have  properties  that 
distinguish  them  from  transition-metal-based  bulk  metallic  glasses  (BMGs),  such  as  very  low 
Young’s  and  shear  moduli  that  are  comparable  with  the  moduli  of  human  bones  [7,8],  very  low 
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density  (1.6 -2. 8  g/cm  )  [9],  and  strong  relaxation  dynamics  of  the  super-cooled  liquid  [10].  Glass 
forming  ability,  as  well  as  physical  and  mechanical  properties  of  metallic  glasses,  is  believed  to 
depend  on  the  type  of  short  range  order  (SRO)  and  medium  range  order  (MRO)  of  alloying 
elements  in  the  amorphous  structure.  It  is  therefore  interesting  to  know  how  the  amorphous 
structure  of  Ca-Mg-based  BMGs  differ  from  transition-metal-based  BMG  structures. 

Structural  analysis  of  metal-metal  BMGs  has  focused  on  transition  metals  glasses  such  as  Zr-Cu 
[  1 1, 12, 13, 14],  Zr-Pt  [15],  and  Zr-Cu-Al  [  16, 17].  The  results  indicate  that  the  packing  of  atoms 
in  these  materials  is  not  random,  but  is  strongly  influenced  by  chemical  interactions.  In 
particular,  icosahedral  SRO  has  been  identified  and  its  importance  on  relaxation  dynamics,  GFA 
and  mechanical  properties  of  these  BMGs  has  been  discussed.  Transition  metal  -  metalloid 
glasses  such  as  Fe-C  and  Ni-B  [  18, 19,20],  which  have  marginal  GFA,  and  many  A1  -  transition 
metal  glasses  with  high  A1  content  [21],  are  not  characterized  by  icosahedral  SRO.  These 
observations  have  been  used  to  support  a  strong  correlation  between  icosahedral  SRO  and  good 
GFA  [22], 

In  spite  of  good  GFA,  the  atomic  structure  of  Ca6oMgxZn4o-x  BMGs  using  X-ray  and  neutron 
diffraction  and  Reverse  Monte  Carlo  (RMC)  simulation  [23]  show  no  icosahedral  SRO.  The 
basic  building  blocks  of  the  amorphous  structure  in  these  glasses  are  Mg-  and  Zn-  centered 
clusters  with  primarily  Ca  atoms  in  the  first  coordination  shell.  The  average  coordination  number 
(CN)  around  Zn  is  9-10  [~6-7  Ca  +  3  (Mg+Zn)],  with  dominant  (0,3, 6,0)  and  (0,2, 8,0) 
polytetrahedral-type  clusters,  and  the  CN  around  Mg  is  1 1-12  [~7-8  Ca  +  4  (Mg+Zn)],  with  a 
large  fraction  of  (0,2,8, 1)  and  (0,2, 8, 2)  clusters.  Ca  atoms  have  an  average  CN  of  13.  Ab  initio 
simulations  support  these  results.  Although  full  icosahedra  do  not  dominate  in  Ca-Mg-Zn  BMGs, 
five-fold  bonds  in  the  form  of  pentagonal  bi-pyramids  are  the  most  populous  structural  units  and 
dominate  in  the  clusters  [23].  A  large  fraction  of  five-fold  bonds  and  the  lack  of  icosahedral  SRO 
has  also  been  found  in  the  molecular  dynamic  (MD)  simulated  amorphous  structure  of  Mg-Cu 
alloys  [24].  The  average  CN  around  Cu  and  Mg  is  composition  dependent,  changing  from  less 
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than  10  with  no  icosahedral  clustering  to  above  12  with  increasing  Cu  content.  Only  on  the  Cu- 
rich  side,  where  the  CN  is  close  to  12,  does  a  predominance  of  five-fold  bonds  indicate  Cu- 
centered  icosaherdal  SRO,  which  is  similar  to  the  Zr-Cu  case  [1 1,12,13,14],  However,  direct 
correlation  between  the  icosahedral  SRO  and  glass  forming  ability  has  not  been  established. 

The  composition  dependence  of  GFA  in  Ca-Mg-Cu  glasses  has  been  well  established  [5,6].  The 
GFA  and  relaxation  dynamics  of  some  of  these  alloys  are  similar  or  even  slightly  better  than 
those  of  the  best  Ca-Mg-Zn  BMGs  [9,10].  At  the  same  time,  Ca-Mg-Cu  BMGs  have  better 
corrosion  resistance  [25,26]  and  seem  to  be  stronger  than  Ca-Mg-Zn  BMGs  [7,8,9].  In  addition, 
strong  covalent-like  interactions  between  Cu-Mg  and  Cu-Ca  have  been  noticed  [5].  It  is  therefore 
useful  to  compare  the  amorphous  structures  of  Ca-Mg-Cu  and  Ca-Mg-Zn  BMGs.  In  the  present 
work,  the  atomic  structures  of  Ca4o+xMg25Cu35-x  BMGs  are  studied  using  a  synergistic 
combination  of  neutron  diffraction,  quantum  molecular  dynamic  (QMD)  simulation  and 
constrained  reverse  Monte  Carlo  (CRMC)  simulation.  These  results  are  compared  with  an  earlier 
analysis  of  the  structure  of  Ca-Mg-Zn  metallic  glasses  [26]. 


2.  EXPERIMENTAL  PROCEDURES 

Fully  amorphous  samples  of  compositions,  Ca4o+xMg25Cu35_x,  where  X  =  0,  5,  10  and  20  at.%, 
were  prepared  by  melt-spinning  in  the  fonn  of  ribbons  as  described  in  [23].  The  densities,  p0,  of 
the  amorphous  alloys  were  measured  with  a  helium  pycnometer  AccuPyc  1330  VI. 03  and  the 
values  (in  g/cm3  and  atoms/A3)  are  given  in  Table  1.  Neutron  diffraction  experiments  were 
conducted  at  room  temperature  under  vacuum  using  the  General  Materials  (GEM)  diffractometer 
at  the  ISIS  high-intensity  pulsed  neutron  source  (Rutherford  Appleton  Laboratory,  Didcot,  UK) 
[27].  GEM  has  eight  detector  hanks  that  collect  data  over  a  wide  range  in  Q  (from  0.1  to  100  A" 

),  where  (9=4;rsin(9/Z  is  the  magnitude  of  the  scattering  vector  for  a  neutron  of  wavelength  X 
scattered  at  an  angle  2  0.  Just  prior  to  the  neutron  diffraction  experiment,  the  ribbon  samples 
were  crushed  into  powder  and  loaded  into  10.3  mm  diameter  cylindrical  vanadium  containers. 
The  thickness  of  the  container  walls  was  25  pm.  Data  reduction,  correction  and  analysis  were 
carried  out  using  the  program  GUDRUN  [28]  and  the  ATLAS  suite  of  programs  [29].  This 
yielded  the  experimental  total-scattering  structure  factors  (SSF),  S(Q)\ 


S(Q )  =  y  ctc AbJAJQ)  - 1]  =  —  —  -  y  C,bf 

7?=\  V  do.  ty 
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Here - is  the  differential  neutron  cross-section  per  unit  solid  angle  Q  for  the  n-  component 

N  dQ 


alloy  (n  =  3  here),  c,  and  b,  are,  respectively,  the  atomic  fraction  and  the  coherent  bound  neutron 
scattering  length  of  element  i,  Ajj(Q)  are  the  Faber-Ziman  partial  structure  factors  [30],  and 

n 

is  the  total-scattering  cross  section  of  the  alloy.  Each  S(Q )  was  Fourier  transformed 

i= 1 

to  give  a  real  space  neutron  total  radial  distribution  function  (RDF),  G(r),  defined  as  [31]: 


n  t  00 

G(r)  =  Z  VM  [?!/  W  - 1]  =  t-2 - \QS(Q)L(Q)sm(rQ)dQ  (2) 

ij= 1  ^  1  Po  0 

Here  pQ  is  the  average  density  of  the  alloy  (in  atoms  per  A3),  r  is  a  distance  from  an  average 
origin  atom  in  the  amorphous  structure,  and  g,j(r)  are  the  partial  radial  distribution  functions 
(PRDF),  which  are  defined  as: 


gy(r)  = 


1  dnv(r) 
Am'2CjP0  dr 


(3) 


where  d//,j  are  the  number  of  elements  of  type  j  between  distances  r  and  r+dr  from  an  element  of 
type  i.  The  Lorch  modification  function  [32],  L(Q),  with  a  maximum  momentum  transfer  Qmax  of 
25  A'1  was  used  to  reduce  tennination  ripples  in  the  Fourier  transform  (Equation  2).  Combining 
equations  (1)  and  (2)  gives  the  following  relation  between  g,j(r)  and  Aj/Q): 

oo 

2*2rp, [«,  ('•)-!]=] QL 4 (6)  -  mQ) sin (rQ)dQ  (4) 

0 

Using  Equation  (3),  the  partial  coordination  number  of  an  element  j  in  the  first  coordination  shell 
of  the  element  i  is  determined: 


Nij  =  4 ncjpo  |  g/~dr  (5) 

Y  • 

'min 

Here  rmin  and  rmax  are  the  positions  of  the  start  and  end  of  the  first  peak  in  the  respective  g,j(r). 
QMD  simulation  of  the  Ca-Mg-Cu  amorphous  structures  was  conducted  using  the  Vienna  Ab- 
initio  Simulation  Package  (VASP)  [33].  The  Projector  Augmented- Wave  (PAW)  method  [34,35] 
and  Perdew-Burke-Emzerhof  (PBE)  functional  [36]  were  used  to  describe  interacting  valence 
electrons.  Cubic  simulation  boxes  with  periodic  boundary  conditions  contained  200  atoms,  and 
the  box  size  was  detennined  by  the  density  of  the  modeled  material.  The  liquid  was  equilibrated 
at  1000K,  then  was  quenched  to  300K  in  100K  temperature  steps,  equilibrating  at  each 
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temperature  for  3000  steps.  Another  3000  consecutive  configurations  were  further  sampled  at 
300K  with  a  time  step  of  5  fs,  and  an  averaged  PRDF  was  calculated  from  these  configurations, 
thus  taking  into  account  thennal  vibrations  and  assuming  no  structural  relaxation  within  15  ps. 
The  QMD  model  provided  simulated  PRDFs,  gijQMD(r),  for  every  studied  alloy. 

To  expand  the  volume  of  the  simulated  structures  and  better  fit  to  experimental  diffraction  data, 
alloy  density  and  closest  approach  constraints,  the  CRMC  simulation  technique  described  in 
detail  in  [37]  was  used.  Simulation  boxes  with  periodic  boundary  conditions  contained  1600 
atoms,  and  the  box  volume  was  detennined  by  the  density  of  the  material.  The  input  data  for  the 
CRMC  simulation  were  S(Q),  G(r),  and  g,yQMD(r).  Renormalization  of  g;/0VII>(r)  was  allowed  to 
correct  possible  QMD  simulation  errors.  More  than  10  million  accepted  random  atom  moves 
occurred  for  each  sample  during  the  CRMC  simulation  before  the  final  configurations  provided 
acceptable  fits  to  the  experimental  data.  This  assured  total  independence  of  the  modeled 
structures  on  the  initial  configurations  of  randomly  distributed  atoms. 

The  simulated  structures  were  statistically  analyzed  using  pair  and  three-body  correlation 
functions,  Voronoi  tessellation  and  nearest  neighbor  approaches,  which  allowed  calculation  of 
local  structural  features  such  as  total  and  partial  CN,  triplet  angle  correlations,  type  and 
distribution  of  characteristic  coordination  polyhedra,  local  packing  fraction  SRO  and  MRO. 

3.  RESULTS 

3,1  Neutron  Diffraction  Analysis 

The  experimental  SSFs,  S(Q),  for  the  four  samples  are  shown  in  Figure  la  and  the  RDFs,  G(r), 
are  shown  in  Figure  2a.  The  S(Q)  curves  damp  very  quickly  with  Q  and  no  oscillations  are 
apparent  beyond  ~15  A"1,  indicating  that  the  glasses  have  a  wide  range  of  interatomic  distances. 
As  the  Cu  concentration  increases,  the  first  sharp  diffraction  peak  (FSDP)  shifts  to  higher  Q  and 
broadens  (Figure  lb).  The  RDFs  show  that  an  increase  in  Cu  content  shortens  the  average 
interatomic  distance,  and  decreases  the  distribution  width  of  the  interatomic  distances  (Figure 
2b).  Thus  the  shift  and  broadening  of  the  FSDP  with  an  increase  in  the  Cu  content  can  be 
interpreted  as  arising  from  the  shortening  of  the  average  interatomic  distance  and  narrowing  of 
the  distribution  of  interatomic  distances,  respectively.  There  is  also  a  pre-peak  at  Q  ~  1.2  A'1  in 
the  S(Q)  curves,  which  indicates  the  presence  of  medium  range  order  in  these  alloys  [38,39]. 
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Systematic  changes  of  the  shape  of  the  first  RDF  peak  with  composition  are  observed  (Figure 
2a).  In  particular,  the  first  RDF  peak  in  Ca6oMg25Cui5  has  two  maxima,  one  at  r  =  3. 1 1  A  and 
another,  higher  intensity,  at  r—  3.63  A.  As  the  amount  of  Cu  increases,  the  relative  intensity  of 
the  first  maximum  increases  and  its  position  slips  to  lower  r-values  (Figure  2b),  while  the  second 
maximum  degenerates  into  a  small  bump  in  alloys  with  25  and  30%  Cu  and  disappears  in  the 
alloy  with  35%  Cu.  At  the  same  time,  an  additional  maximum  starts  to  develop  at  r  ~  2.63  A, 
preceding  the  first  maximum,  in  the  alloys  with  30%  and  35%  Cu.  Such  a  strong  concentration 
dependence  of  the  shape  of  the  first  RDF  peak  is  apparently  due  to  an  increase  in  the  fraction  of 
shorter  Cu-Cu  and  Cu-Mg  bond  pairs  and  a  decrease  in  the  fraction  of  longer  Ca-Ca  and  Ca-Mg 
bond  pairs  with  an  increase  in  the  amount  of  Cu. 

3.2  Quantum  Molecular  Dynamic  Simulation 

Figure  3  shows  the  six  QMD-simulated  PRDFs  for  the  Ca40Mg25Cu35,  Ca50Mg25Cu25  and 
Ca6oMg25Cui5  alloys.  The  PRDFs  are  rather  smooth  since  they  account  for  atom  vibrations.  The 
effect  of  composition  on  the  width  and  position  of  the  first  PRDF  peak  is  weak  relative  to  the 
effect  on  the  positions  and  intensities  of  the  higher  order  peaks.  This  indicates,  to  the  first  order, 
that  bond  distances  between  respective  atom  pairs  in  the  first  coordination  shell  are  not  affected 
by  composition  and  that  the  partial  coordination  number  of  element  j  around  element  i  is  linearly 
proportional  to  the  concentration  Cj  and  alloy  density  pQ  (see  Equations  3  and  5).  From  the  QMD- 
simulated  PRDFs  and  using  Equation  (2),  simulated  total  RDFs  were  obtained  for  every  alloy 
and  compared  with  the  respective  experimental  neutron  RDFs  (Figure  4).  A  very  good  match  is 
seen.  Since  QMD  simulations  were  done  completely  independently  from  the  neutron  diffraction 
experiments,  these  good  fits  indicate  that  the  simulated  PRDFs  are  acceptable  and  reliable. 

The  average  partial  and  total  coordination  numbers  around  Ca,  Mg  and  Cu  atoms,  from  the  QMD 
simulations,  are  given  in  Table  2  and  the  most  frequent  (mode,  r™ode ),  weighted  average  (mean, 

f.Mea"  j  anc|  cut-off  (maximum,  rf"'-""  )  distances  between  atom  pairs  in  the  first  coordination 

shell  are  given  in  Table  3,  Table  4,  and  Table  5,  respectively.  The  weighted  average  distances, 
yMec," ,  between  the  i  and  j  atoms  were  calculated  using  the  following  equation  [23]: 

rCut-off  I  rCut-off 

r™ean  =  |  rijgij{r)dr  \gij(r)dr  (6) 

0  /  0 
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3.3  Constrained  Reverse  Monte  Carlo  Simulation 

The  QMD-simulated  PRDFs,  g^MD(r),  were  used  as  soft  constraints  in  CRMC  simulations  to 
force  the  local  atomic  order  to  be  consistent  with  experimental  S(Q)  and  G(r)  and  simulated 
g,/QMD(r).  Renormalization  of  g,/QMI)(r)  during  CRMC  allowed  minor  adjustment  of  the  QMD 
simulation  to  give  much  better  fits  of  the  CRMC-simulated  S(Q)  and  G(r )  to  the  experimental 
functions  (Figure  5).  The  CRMC  PRDFs  are  shown  in  Figure  6.  Similar  to  QMD  results,  the  first 
peak  position  almost  does  not  depend  on  alloy  composition,  however,  the  peak  intensity  does. 
Using  these  PRDFs  and  Equation  (3),  the  partial  coordination  numbers,  as  well  as  the  mode, 
mean  and  maximum  cut-off  distances  between  the  pair  atoms  in  the  first  coordination  shell,  were 
calculated  and  the  results  are  given  in  Table  2,  Table  3,  Table  4,  and  Table  5,  respectively. 

3.3.2  Voronoi  Analysis  of  Ca-Mg-Cu  Amorphous  Structure 

Voronoi  analysis  was  used  to  study  local  atom  environments  beyond  the  PRDFs  results  [23,40]. 
In  Voronoi  analysis,  atoms  that  have  common  faces  in  their  Voronoi  polyhedra  are  nearest 
neighbors  and  a  coordination  polyhedron  with  vertices  at  the  first-shell  atom  positions  and  edges 
coinciding  with  the  interatomic  bonds  is  defined  for  any  chosen  atom  /  [41].  Each  coordination 
polyhedron,  also  called  an  /-centered  cluster  [40,42,43],  is  associated  with  the  respective  /- 
centered  Voronoi  polyhedron  and  can  be  assigned  a  Voronoi  signature  (v3,v4,v5,v6).  For  Voronoi 

polyhedra,  vm  is  the  number  of  faces  containing  m  edges;  while  for  the  respective  coordination 

polyhedron,  vm  is  the  number  of  vertices  common  to  m  polyhedron  edges  (or  faces)  [44].  In  the 

latter  case,  m  is  also  called  the  vertex  coordination.  Those  /-centered  coordination  polyhedra  with 
the  same  Voronoi  signature  are  considered  to  be  topologically  equivalent  (even  though  they  may 
not  be  identical)  because  they  can  be  transformed  into  each  other  without  changing  the  number 
of  vertices  and  connecting  edges.  In  addition  to  assigning  topology,  the  Voronoi  signature  also 
defines  the  total  coordination  number  (/V,j  of  the  /-centered  cluster  as  N,  =  Y.vm.  Clusters  with  the 

same  Voronoi  signature  can  be  chemically  different.  Therefore,  in  addition  to  the  Voronoi 
signature,  partial  coordination  numbers,  i.e.  the  number  of  atoms  of  different  species,  should  also 
be  known  for  a  more  complete  description  of  the  SRO  in  the  amorphous  structure.  The  partial 
and  total  coordination  numbers  obtained  with  Voronoi  tessellation  of  the  CRMC-simulated 
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amorphous  structure  of  Ca4o+xMg25Cu35_x  alloys  are  given  in  Table  2.  They  are  in  good 
agreement  with  those  obtained  with  the  use  of  the  CRMC-simulated  PRDFs. 

Figure  7  illustrates  the  types  and  fractions  of  Ca-,  Mg-  and  Cu-  centered  coordination  polyhedra 
in  the  CRMC-simulated  amorphous  structures.  Although  many  types  of  coordination  polyhedra 
are  present,  the  most  common  clusters  are  (0,2, 8, 4)  for  Ca-centered,  (0,2, 8, 2)  for  Mg-centered, 
and  (0,3, 6,0)  and  (0,2,8, 1)  for  Cu-centered  clusters.  The  fractions  and  distributions  of  the  clusters 
depend  on  alloy  composition.  Figure  8  shows  the  fractions  of  3-,  4-,  5-,  and  6-  coordinated 
vertices  in  the  Ca-,  Mg-  and  Cu-  centered  clusters  in  the  Ca-Mg-Cu  amorphous  alloys.  The  five- 
coordinated  vertices  dominate  in  all  clusters,  which  is  an  indication  that  pentagonal  bi-pyramids 
are  the  main  building  blocks  in  these  amorphous  alloys.  Taking  into  account  that  the  pentagonal 
bi-pyramid  consists  of  five  tetrahedra,  one  may  conclude  that  polytetrahedral  packing  plays  an 
important  role  in  the  formation  of  the  amorphous  structure  in  the  Ca-Mg-Cu  alloys. 

Radical  Voronoi  tessellation  gives  the  Voronoi  cell  volume  (VCV)  of  each  individual  atom.  The 
average  VCV  values  of  Ca,  Mg  and  Cu  atoms  are  given  in  Table  6.  The  VCV  is  higher  for  larger 
atoms.  The  VCV  of  Ca  and  Cu  increases,  while  that  of  Mg  is  nearly  constant  or  increases 
slightly,  with  decreasing  Cu  concentration.  Assigning  the  effective  atomic  radius  of  element  i  as 
half  the  mode  i-i  bond  distance  (Table  3),  the  average  void  volume  and  packing  fraction  near  i 
atoms  were  calculated  for  each  alloy  as  the  difference  between  the  respective  Voronoi  volume 
and  the  volume  of  atom  /,  and  as  the  fraction  of  the  atom  volume  relative  the  VCV,  respectively 
(Table  6).  The  volume  of  voids  and  the  packing  fraction  are  higher  near  larger  atoms.  The  void 
volume  increases  but  the  packing  fraction  decreases  near  Ca  and  Mg  atoms  with  a  decrease  in  Cu 
concentration.  The  packing  fraction  of  Cu  does  not  depend  on  the  alloy  composition  within  the 
25-35%  Cu  range.  However,  it  decreases  in  Ca6oMg25Cui5.  The  average  packing  fraction  of  the 
alloys  determined  from  their  density  (Table  6)  shows  similar  trends  and  increases  from  0.68  to 
0.73  when  the  Cu  content  increases  from  15%  to  35%.  These  values  show  that  packing  is 
significantly  more  efficient  than  in  monatomic,  dense-random  packed  structures,  and  is 
comparable  to  packing  efficiencies  in  crystalline  metals. 

3.3.3  Three-Body  Correlations 

Three-body  correlations  can  provide  additional  information  about  the  most  common  local 
structures  around  a  centered  atom  because  the  characteristic  angles  between  bonds  connecting 
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any  two  atoms  located  in  the  first  or  second  coordination  shell  with  the  centered  atom  are  very 
sensitive  to  the  spatial  atom  distribution  [22].  Figure  9  illustrates  the  bond  angle  distributions  in 
the  Ca45Mg25Cu3o  amorphous  alloy  within  the  first  and  second  coordination  shells.  The  upper 
limits  of  the  bond  lengths  for  the  first  shell  were  set  to  the  cut-off  values  given  in  Table  5,  while 
the  bond  length  ranges  for  the  second  shell  were  set  between  the  values  corresponding  to  the  first 
and  second  minima  of  the  PRDFs  (see  Figure  6).  Two  characteristic  maxima,  one  in  the  range  of 
50-70°  and  another  between  100-140°,  are  clearly  seen  in  the  triplet  distributions  for  the  first 
coordination  shell  (Table  7).  At  least  three  maxima  were  identified  in  bond  angle  distributions 
for  the  second  coordination  shell.  These  are  near  30°,  60°  and  90°  (Table  8). 

3.  DISCUSSION 

QMD  simulated  structures  give  smaller  CNs  for  Cu-Cu,  Cu-Mg,  Mg-Mg  and  Mg-Cu  and  larger 
Ca-Ca,  Ca-Mg,  Ca-Cu,  Mg-Ca  and  Cu-Ca  CNs  than  does  CRMC  (Table  2).  QMD  simulated 
structures  also  give  larger  mean  (Table  4)  and  cut-off  (Table  5)  pair  atom  distances.  Given  the 
very  high  QMD  quench  rate  and  the  periodic  boundary  condition  on  a  relatively  small  volume, 
these  minor  discrepancies  in  structural  details  are  not  surprising.  At  the  same  time,  CNs  obtained 
with  Voronoi  tessellation  are  in  full  agreement  with  CNs  obtained  directly  from  the  CRMC 
partial  RDF’s.  Thus,  only  the  CRMC-simulated  amorphous  structures  will  be  discussed  here. 

3.1.  Pair  Bond  Distance  Correlations 

fMode  values  (Table  3)  are  always  smaller  than  the  respective  rMea„  values  (Table  4).  The 
difference  is  smallest  for  Ca-Ca  pairs  (3.2%)  and  increases  in  the  order  of  Ca-Mg  (4.3%),  Mg- 
Mg  (7.7%),  Ca-Cu  (8.2%),  Cu-Cu  (8.8%),  and  Mg-Cu  (13.5%).  This  indicates  that  the  first-shell 
gij(r)  peaks  are  asymmetric  (non-Gaussian),  especially  for  Cu-  containing  pairs.  rMode  and  rMean 
values  can  be  compared  with  atomic  separations  in  crystalline  metals  and  alloys.  Crystalline 
metallic  separations,  rM,  and  covalent  distances,  r(--,  were  estimated  as  a  sum  of  metallic  and 
covalent  radii  of  respective  elements  (Table  9).  The  metallic  radii  are  half  the  nearest  distance 
between  atoms  in  crystal  lattices  of  the  respective  pure  metals  [69].  The  covalent  atomic 
separations  [70]  were  obtained  from  analysis  of  nearest  distances  in  crystalline  intennetallic 
compounds  using  the  Cambridge  Structural  Database.  From  Table  9,  the  rMea,i  and  rMode  values 
for  Ca-Ca  atom  pairs  are  smaller  than  rM,  so  that  the  mean  and  mode  Ca-Ca  separations  in  Ca- 
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Mg-Cu  glasses  are  slightly  compressed.  However,  for  all  other  atomic  pairs  rMean  values  are 
larger  and  rMode  values  are  smaller  than  rM.  The  larger  mean  separations  suggest  that  some  Mg 
and,  especially,  Cu  atoms  become  slightly  displaced  from  each  other.  It  is  likely  that,  similar  to 
Zn  in  Ca-Mg-Zn  amorphous  alloys  [23],  some  Cu  atoms  occupy  gaps  between  the  larger  Ca  and 
Mg  atoms  in  the  first  coordination  shell  and  their  distances  from  the  center  atom  vary  depending 
on  the  gap  sizes.  At  the  same  time,  rMode  for  the  atom  pairs  containing  Cu  are  even  smaller  than 
the  respective  covalent  bond  distances  in  crystalline  compounds,  and  the  mode  separations  of 
Ca-Cu  and  Mg-Cu  are  well  represented  by  the  comparison  between  rMode  and  re-  These  data 
show  that  the  absence  of  long-range  order  constraints  allows  shortening  of  first  neighbor 
interatomic  distances  in  amorphous  structures,  as  compared  to  bond  distances  in  crystalline 
metals  and  alloys. 

The  Mg-Cu  (2.67  A)  and  Ca-Cu  (3.04  A)  mode  pair  distances  and  are  much  shorter  than  those 
calculated  from  the  Ca-Ca,  Mg-Mg  and  Cu-Cu  rMode  values  (i.e.  2.81  A  and  3.13  A,  respectively, 
see  Table  9).  Shortening  of  Mg-Cu  and  Y-Cu  bonds  has  been  observed  in  amorphous 
MgeoCusoYio  [45].  Bond  shortening  has  been  discussed  in  amorphous  and  crystalline  systems 
[46,47,48],  and  is  likely  related  to  sp-d  electron  hybridization  [49].  This  shortening  is  often 
associated  with  a  covalent  bonding  component  and  indicates  strong  chemical  interactions, 
suggesting  that  attractive  forces  will  increase  in  the  sequence:  Ca-Mg,  Ca-Cu  and  Mg-Cu.  It  has 
been  suggested  that  the  thermodynamic  heat  of  mixing  can  represent  this  chemical  interaction 
[50],  and  experimental  data  [51,52,53]  give  the  heats  of  mixing  of  equimolar  binary  alloys  as 
-6.0±1.0  kJ/mol  for  Ca-Mg,  -6.5±1.0  kJ/mol  for  Ca-Cu  and  -9.0±1.0  kJ/mol  for  Cu-Mg.  These 
values  give  the  same  trend  as  for  bond  shortening,  but  the  link  between  heats  of  mixing  and  bond 
lengths  is  phenomenological.  Further,  the  heats  of  mixing  give  a  global  representation  of 
interatomic  bond  enthalpies,  which  include  the  number  of  bonds  formed  and  the  energies  of 
bonds  between  two  atoms.  Both  of  these  values  are  likely  to  depend  upon  composition  and  local 
atom  environment,  so  that  trends  in  relative  bond  lengths  obtained  from  heats  of  mixing  should 
be  used  with  caution. 

The  observation  that  nearest-neighbor  bond  lengths  are  shorter  than  those  in  the  competing  long- 
range  ordered  crystals  may  help  explain  why  these  are  easy  glass  fonners.  The  equilibrium 
crystalline  phases  that  compete  with  the  amorphous  phase  appear  to  minimize  total  energy  by 
sacrificing  optimal  short-range  order  to  achieve  beneficial  long-range  order.  On  the  other  hand, 
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the  absence  of  long-range  atomic  order  allows  metallic  glasses  to  reduce  the  energy  difference 
between  the  metastable  glass  and  equilibrium  crystal(s)  by  optimizing  short-range  atomic 
interactions  and  arranging  the  atoms  in  efficiently  packed  clusters  [40,42,54,55].  The  largest 
contribution  to  condensed  phase  stability  is  expected  to  come  from  nearest-neighbor  interactions 
and  volume  minimization,  so  that  the  optimized  short-range  interactions,  as  well  as  the  presence 
of  MRO,  in  metallic  glasses  can  produce  structures  that  have  only  a  small  energetic  disadvantage 
relative  to  the  equilibrium  crystalline  structure  [22,23].  Kinetic  constraints  from  quenching 
restrict  the  long  range  atomic  redistribution  needed  to  achieve  long-range  order  and  to  further 
minimize  the  total  system  energy,  thus  favoring  glass  fonnation.  The  present  results  suggest  that 
metallic  glasses  not  only  have  short  range  order,  but  may  have  “better”  combination  of  short- 
range  topological  and  chemical  order  (in  tenns  of  optimal  bond  length  and/or  atomic 
arrangement)  than  the  competing  crystals. 

3.2  Coordination  Numbers  and  Local  Chemical  Order 

Increasing  the  Cu  concentration  from  15  to  35  at%  continuously  increases  the  total  coordination 
number  around  Ca,  CNca,  from  13.6  to  15.0  due  to  a  faster  increase  in  the  number  of  smaller  Cu 
atoms  ( CNca-Cu  increases  from  1.7  to  4.4)  and  a  slower  decrease  in  the  amount  of  the  larger  Ca 
atoms  (CN Ca-Ca  decreases  from  8.6  to  6.9)  (Figure  1 1).  On  the  other  hand,  the  CNMg  and  GVq, 
almost  do  not  depend  on  alloy  composition.  However,  similar  to  the  environment  around  Ca, 
CNms-Cu  and  CNcu-Cu  increase  from  1.4  to  3.7  and  from  0.6  to  2.6,  respectively,  while  CNMg-Ca 
and  CNcu-Ca  decrease  from  8.0  to  5.9  and  from  6.8  to  5.0,  respectively,  with  the  increase  in  Cu 
concentration.  The  number  of  Mg  atoms  around  Ca,  Mg  and  Cu  very  weakly  depends  on  alloy 
composition  and,  on  average,  CNca-Mg  =  3.5±0.2,  NMg-Mg  =  2.8±0.3  and  CNcu-Mg  =  2.7±0.3  (see 
Table  2  and  Figure  11),  which  may  indicate  a  rather  homogeneous  Mg  distribution.  The 
dependences  of  the  partial  coordination  numbers  on  Cu  concentration  can  be  described  by  linear 
equations  (Table  10).  Extrapolation  to  low  Cu-concentrations  predict  no  Cu-Cu  first-neighbor 
interactions  in  alloys  with  Cu  concentrations  less  than  10.3  at.%.  At  these  low  Cu  concentrations, 
CNcu-Ca  and  CNcu-Mg  are  predicted  to  be  7.3  and  2.5,  respectively,  CNca  =  13  and  CVm<,  =  1 1.8. 

To  detennine  the  degree  of  chemical  short  range  order  (CSRO)  around  Ca,  Mg  and  Cu  atoms, 
the  Warren-Cowley  parameter,  a,,,  which  is  defined  as  [56]: 
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ccij  =  1 CNij/(cjCNi)  (7) 

was  used.  Here  Cj  is  the  atomic  fraction  of  the  element  j,  CNj  is  the  partial  coordination  number 
of  element  j  around  the  element  i,  and  CN,  is  the  total  coordination  number  around  element  i.  A 
negative  ay  (i  ±  j )  indicates  the  presence  of  CSRO  (i.e.  the  number  of  atoms  j  in  the  first  shell 
exceeds  the  average  concentration).  For  a  random  solution,  ay  are  zero.  A  modified  CSRO 
parameter,  ayjk),  which  is  defined  as 

ocifjk)  =  1  -  (CNij+  CN ,;k)/ {{cj+Ck)CN i)  for  i  +  k  (8) 

can  also  be  used  to  analyze  the  CSRO  in  ternary  alloys.  Negative  a^k)  values  indicate  the 
presence  of  CSRO,  while  positive  values  indicate  the  presence  of  chemical  short  range  clustering 
(CSRC)  near  the  i  atom,  a^k)  and  an  are  dependent  parameters,  namely 

(1-C;)  CXi(jk)  Cj  Ctii  (9) 

Therefore,  a„  <  0  should  indicate  CSRC  (i.e.  increased  local  concentration  of  like  atoms)  and  a„ 
>  0  should  indicate  CSRO  around  atom  i. 

Values  of  a,,  and  ai(jk)  for  Ca-Mg-Cu  glasses  are  given  in  Table  1 1  and  Table  12,  respectively. 
The  alloys  show  CSRC  of  Ca  atoms  that  is  largest  for  the  alloy  with  40%  Ca  (acaCa  =  -0.16)  and 
it  tends  to  decrease  with  an  increase  in  Ca  concentration  (at  60%  Ca  acaCa  =  -0.05).  Pronounced 
CSRO  is  detected  near  Cu  atoms.  These  local  regions  are  enriched  with  Ca  and  Mg,  supporting 
the  earlier  observation  of  strong  Ca-Cu  and  Mg-Cu  bonding.  Near  Mg  atoms,  oiMgCa  is  negative 
and  aMgcu  is  positive,  indicating  enrichment  in  Ca  and  a  Cu  deficit  in  Mg-centered  clusters.  The 
a.ijk)  criterion  shows  clear  evidence  of  CSRC  near  Ca  atoms  and  CSRO  near  Cu  atoms  when 
unlike  elements  are  considered  together  (Table  12).  At  the  same  time,  a^k)  is  almost  zero  or 
slightly  negative  for  Mg-centered  clusters  indicating  neutral  environment. 


3.3.  Voronoi  Tessellation  Analysis 

Voronoi  analysis  indicates  that  Kasper- type  polyhedra  dominate  Ca-Mg-Cu  glass  structures. 
Among  the  clusters  in  Figure  7,  (0,4, 4,0),  (0,3, 6,0),  (0,2, 8,0),  (0,2,8, 1),  (0,0,12,0),  and  (0,1,10,2) 
are  non-distorted  Kasper  polyhedra,  whereas  (0,3,6, 1),  (0,4, 4, 2),  (0,3, 6, 2),  (0,4, 4, 3),  (0,2, 8, 2), 
(0,3, 6, 3),  (0,4, 4, 4),  (0,2, 8, 3),  (0,3, 6, 4),  (0,4, 4, 5),  (0,1,10,3),  (0,2, 8, 4),  (0,3, 6, 5),  (0,4, 4, 6), 


(0,1,10,4),  (0,2, 8, 5),  (0,3, 6, 6),  are  distorted  Kasper  polyhedra  with  four-fold  and  six-fold 
disclinations  [22],  and  only  (0,3, 7, 4),  (1,2, 6, 3),  (1,2, 5, 2),  (1,3, 3, 2)  and  (1,2, 5, 3)  are  other-type 
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polyhedra.  The  Kasper  polyhedra  edges  can  have  different  lengths-  this  is  inevitable  in  glass; 
therefore,  they  are  not  identical  even  though  they  share  the  same  topology  (in  terms  of  the 
Voronoi  index)  with  the  latter.  Similar  to  other  work  [22],  a  cluster  is  identified  here  as  a  non- 
distorted  Kasper  polyhedron  if  it  has  a  Voronoi  index  of  the  Kasper  polyhedron  and  the  tenn 
‘distorted’  is  assigned  to  Kasper  polyhedra  with  disclinations. 

Kasper  polyhedra  and  their  distorted  variants  account  for  over  50%  of  the  nearest-neighbor 
clusters  (Table  13).  Since  Kasper  polyhedra  are  polytetrahedral,  and  many  other  polyhedra  also 
contain  tetrahedra,  we  conclude  that  the  topological  SRO  of  Ca-Mg-Cu  is  polytetrahedral  in 
nature,  which  has  been  demonstrated  in  many  other  metallic  glasses  [22].  The  current  study 
therefore  further  supports  the  idea  that  the  SRO  of  MGs  is  characterized  by  polytetrahedral 
packing  via  Kasper  clusters  and  their  distorted  variants.  The  regularity  of  polytetrahedral  packing 
can  be  evaluated  using  the  ratio  between  the  number  of  regular  Kasper  polyhedra  and  total 
number  of  regular  and  distorted  Kasper  polyhedra.  We  find  that  the  atomic  packing  around 
smaller  atoms  (Cu)  is  much  more  regular  polytetrahedral  than  around  larger  atoms  (Ca).  For 
example,  the  fractions  of  regular  Kasper  polyhedra  around  Cu  is  about  55%,  around  Mg  range 
from  23%  to  32%  and  decreases  to  13-20%  for  Ca-centered  clusters.  This  is  reminiscent  of  Cu- 
Zr,  where  the  packing  around  Cu  appears  much  more  regular  than  the  packing  around  Zr  [22], 
The  more  regular  packing  around  Cu  can  indicate  that  Cu-centered  clusters  are  the  primary 
structure-forming  clusters,  whereas  Mg-  and  Ca-centered  clusters  are  secondary  clusters 
resulting  from  specific  arrangements  of  atoms  in  the  vertices  of  and  voids  between  the  Cu 
clusters  (see  section  3.5  below).  This  can  also  be  explained  by  very  strong  interactions  of  Cu 
with  Ca  and  Mg  resulting  in  bond  shortening.  This  feature  could  be  quite  general  for  metallic 
glasses,  i.e.,  polytetrahedral  packing  is  more  prominent  and  better  established  around  smaller 
atoms.  This  is  probably  because  smaller  atoms  have  smaller  CNs,  and  thus  have  limited  choices 
for  efficient  packing,  in  comparison  to  larger  atoms  with  more  neighbors  and  possibilities. 
Similar  to  Ca-Mg-Zn  BMGs  [23],  the  amorphous  structures  of  Ca-Mg-Cu  glasses  contain  a  very 
low  fraction  of  icosahedral  SRO,  in  spite  of  the  observation  that  five-coordinated  vertices 
dominate  in  all  clusters.  A  large  fraction  of  five-fold  bonds  and  the  lack  of  icosahedral  SRO  have 
also  been  found  in  the  MD-simulated  amorphous  structure  of  Mg-Cu  alloys  [24],  The  dominance 
of  five-fold  vertices  suggests  that  tetrahedra  prefer  to  cluster  into  pentagonal  bi-pyramids. 
Analysis  of  amorphous  Mg6oCu3oY io  [Error!  Bookmark  not  defined.]  shows  a  prevalence  of 
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fivefold  bonds,  with  local  atom  arrangements  very  similar  to  the  competing  Mg2Cu  (Laves 
phase)  and  YCm  crystals.  Similar  comparison  of  coordination  polyhedra  in  the  amorphous 
structure  of  Ca6oMgxZn4o-x  with  polyhedra  in  the  competing  crystal  phases,  CaMgi  and  CaZn 
[57],  has  shown  that  SRO  arrangements  around  Ca  and  Mg  atoms  were  different  from  those  in 
the  crystal  phases  [23],  However,  the  (0,3, 6,0)  Zn-centered  coordination  polyhedron,  typical  to 
CaZn,  was  one  of  the  most  common  in  the  amorphous  structure. 

Amorphization  of  CaeoMg^sCuis  competes  with  CaMg?  compound  formation,  amorphization  of 
Ca4oMg25Cu35  competes  with  the  Cu2Mg  phase,  while  solidification  of  Ca45Mg25Cu3o  and 
Ca50Mg25Cu25  alloys  forms  three  equally  probable  crystal  phases,  CaMg2  CmMg  and  CaCu 
[5,6].  CaMg2  is  a  hexagonal  Laves  phase  (space  group  P63/mmc,  space  group  number  194)  with 
lattice  parameters  a  =  5.170  A  and  c  =  8.50  A.  This  crystal  structure  contains  three  coordination 
polyhedra,  two  of  which  are  Mg-centered  (0,0,12,0)  icosahedra  with  6  Mg  and  6  Ca  atoms  in  the 
first  shell,  and  the  third  is  a  Ca-centered  (0,0,12,4)  Friauf  polyhedron  containing  12  Mg  and  4  Ca 
atoms  in  the  first  shell.  CmMg  is  a  cubic  Laves  phase  (space  group  is  Fd-3m,  space  group 
number  227)  with  the  lattice  parameter  a  =  6.990  A.  The  crystal  structure  of  this  phase  has  two 
coordination  polyhedra:  a  Cu-centered  (0,0,12,0)  icosahedron  with  6  Cu  and  6  Mg  atoms  in  the 
first  shell,  and  a  Mg-centered  (0,0,12,4)  Friauf  polyhedron  with  12  Cu  and  4  Mg  atoms  in  the 
first  shell.  The  CaCu  phase  has  a  primitive  crystal  structure  (space  group  P2i/m,  space  group 
number  1 1)  with  lattice  parameters  a  =  19.47  A,  b  =  4.27 1  A,  c  =  5.880  A,  and  p  =  94.30°.  It  has 
10  characteristic  coordination  polyhedra,  5  of  which  are  Cu-centered  Kasper  polyhedra  (0,3, 6,0) 
with  2  Cu  and  7  Ca  atoms  in  the  first  shell,  and  the  other  5  are  Ca-centered  polyhedra  (3, 1,6, 6,1), 
(6, 1,2, 2, 4, 0,2),  (1,3, 6, 4, 3),  (3, 3, 3, 4, 3, 0,1)  and  (1,3, 7, 6),  with  7  Cu  and  10  Ca  atoms  in  the  first 
shell.  [58]  (Additional  indices  in  the  Voronoi  signatures  of  the  Ca-centered  polyhedra  given 
above  correspond  to  the  number  of  7-,  8-,  and  9-  coordinated  vertices,  respectively.) 

None  of  the  Ca-centered  coordination  polyhedra  of  the  competing  crystal  phases  are  present  in 
amorphous  structures  of  Ca-Mg-Cu.  The  icosahedral  SRO  in  the  CaMg2  and  Cu2Mg  crystal 
phases  near  Mg  and  Cu  atoms,  respectively,  is  not  typical  of  Ca-Mg-Cu  glass  structures.  The 
(0,3, 6,0)  Cu-centered  coordination  polyhedron  in  the  CaCu  phase  is  also  common  in  the  studied 
amorphous  structures.  This  (0,3, 6,0)  Cu-centered  cluster  is  probably  the  only  link  between  the 
competing  crystal  and  amorphous  states  and,  therefore,  a  reduced  fraction  of  this  cluster  in  the 
amorphous  structures  should  indicate  better  GFA.  Indeed,  the  Ca6oMg25Cui5  and  Ca40Mg25Cu35 
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amorphous  alloys  have  much  higher  fraction  of  the  (0,3, 6,0)  Cu-centered  clusters  than  the  better 
glass  fonning  Ca45Mg25Cu3o  and  Ca50Mg25Cu25  alloys. 

3.4.  Three-Body  Correlations 

The  Voronoi  analysis  above  indicates  that  tetrahedra  are  dominant  building  blocks  in  amorphous 
structures  of  Ca-Mg-Cu  alloys.  With  more  than  50%  frequency,  the  tetrahedra  form  pentagonal 
bipyramids  and  Kasper  polyhedra.  Because  the  characteristic  vertex  angles  for  regular  tetrahedra 
are  60°  and  for  the  regular  pentagonal  bipyramids  are  60°  and  108°,  the  deviations  of  angle  peak 
positions  in  three-body  correlations  indicate  distortions  from  the  regular  configurations.  These 
distortions  are  mainly  caused  by  different  radii  of  Ca,  Mg  and  Cu  atoms,  as  well  as  non-Gaussian 
distributions  of  the  bond  distances  in  the  first  coordination  shell. 

Table  7  gives  the  first  and  second  maxima  in  bond  angle  distributions  of  Ca-Ca-Ca,  Ca-Mg-Ca, 
Ca-Cu-Ca,  Mg-Mg-Mg  and  Cu-Cu-Cu  triplets  in  the  first  coordination  shell.  If  like  atoms  are 
packed  closely,  then  the  characteristic  angle  for  the  i-i-i  triplets  will  be  equal  or  close  to  60°.  The 
first  peak  for  Ca-Ca-Ca,  Mg-Mg-Mg  and  Cu-Cu-Cu  triplets  occurs  at  -57.9°,  56.8°  and  56.2°, 
respectively.  The  peak  shift  to  angles  less  than  60°  indicates  that  the  i-i  distances  within  the  first 
shell  are  longer  than  the  bonds  between  the  center  and  two  first-shell  i  atoms.  From  Section  3.1, 
we  take  the  distances  between  the  center  and  shell  like  atoms  as  rMode  and  the  distance  between 
the  shell  atom  pair  as  rMean,  so  that  the  angle  of  an  i-i-i  triplet  is  estimated  as 

a  =  COS“‘(r Meaner Mode)  (10) 

The  calculated  angles  for  Ca-,  Mg-  and  Cu-triplets  using  averaged  rMode  and  rMean  values  from 
CRMC  simulations  (see  Table  9)  are  58.9°,  57.4°  and  57.0°,  respectively.  These  are  all  within  1° 
of  measured  angles,  giving  good  agreement. 

For  Ca-Mg-Ca  and  Ca-Cu-Ca  triplets,  the  first  peaks  are  at  ~62.4°  and  69.4°,  respectively.  These 
larger  angles  result  from  the  smaller  atomic  radii  of  Mg  and  Cu  atoms.  Using  rMode  for  Ca-Ca, 
Ca-Mg  and  Ca-Cu  separations,  the  Ca-Mg-Ca  and  Ca-Cu-Ca  triplet  angles  are  estimated  as  63.0° 
and  67.5°,  respectively.  This  is  very  good  agreement  for  Ca-Mg-Ca,  and  shows  reasonable 
agreement  for  Ca-Cu-Ca  given  the  uncertainty  in  bond  angles. 

The  second  bond-angle  distribution  peak  for  Ca-Ca-Ca  triplets  is  located  near  101°  in  the 
Ca4oMg25Cu35  alloy  and  shifts  to  -106°  in  Ca60Mg25Cui5.  This  is  a  little  smaller  than  the 
108°interior  angle  of  a  regular  pentagon,  and  may  be  due  to  distortion  when  some  vertices  are 
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occupied  by  smaller  Mg  and  Cu  atoms.  This  also  explains  the  angle  decrease  with  increasing  Cu 
concentration.  Due  to  shorter  Ca-Cu  and  Ca-Mg  bonds,  the  interior  angle  at  the  Cu  or  Mg  vertex 
increases  at  the  expense  of  other  interior  angles,  the  sum  of  which  is  constant  at  540°.  Indeed,  the 
second  peak  for  Ca-Mg-Ca  and  Ca-Cu-Ca  triplets  is  located  at  1 1 1-1 14°  and  at  121-130°, 
respectively.  Using  these  values  of  interior  angles  at  Ca,  Mg  and  Cu  vertices,  on  average  2.6  Ca, 
1.1  Mg  and  1.3  Cu  atoms  in  the  Ca40Mg25Cu35  alloy  and  ~4.0  Ca,  0.7  Mg  and  0.3  Cu  atoms  in 
the  Ca6oMg25Cui5  alloy  are  estimated  per  pentagon  to  satisfy  the  sum  of  vertex  angles  of  540°. 
The  first  three  maxima  in  bond  angle  distributions  of  Ca-Ca-Ca,  Ca-Mg-Ca,  Ca-Cu-Ca,  Mg-Mg- 
Mg  and  Cu-Cu-Cu  triplets  in  the  second  coordination  shell  are  given  in  Table  8  and  Figure  9b.  If 
like  atoms  contact  in  the  second  shell,  then  the  characteristic  first  angle  for  i-i-i  triplets  will  be 
equal  or  close  to  30°.  This  angle  is  estimated  for  Ca-,  Mg-  and  Cu  triplets  by  assuming  the  i-i 
pair  in  the  2nd  coordination  shell  are  separated  by  rmode  and  that  the  second  shell  is  displaced  from 
the  central  i  atom  by  a  distance  (rMode  +  r),  where  r  is  the  weighted  average  of  rmode  values  for 

the  atoms  in  the  first  shell  of  the  i  atom.  The  triplet  angle  is  estimated  as 


P  c°s  1  |l-(l  2 )[>;*„*.  ('W*. + f 


)!} 


(ii) 


The  angles  for  a  given  i-i-i  triplet  are  averaged  over  all  four  compositions.  The  estimated  Ca-, 
Mg-  and  Cu  triplet  angles  are  30.8°,  27.9°  and  24.4°,  respectively,  and  are  all  within  1°  of  the 
measured  angles  (30.3°,  27.3°  and  24.5°).  A  similar  topological  analysis  of  Ca-Mg-Ca  and  Ca- 
Cu-Ca  triplets  for  the  2nd  coordination  shell  give  estimated  angles  of  32.2°  and  34.0°,  which 
compare  well  with  the  measured  values  of  3 1 .5°  and  32.8°. 

A  second  peak  is  expected  near  60°  when  a  full  atom  occurs  between  the  i-i  pair  in  the  second 
shell,  and  the  bond  angle  distributions  show  a  rather  broad  peak  centered  very  near  60°  for  the 
five  triplets  studied  (Figure  9b).  To  model  this,  an  atom  with  the  effective  radius  of  atoms  in  the 


1SI  shell  of  i  atoms  is  placed  between  the  i-i  atom  pair  in  the  2nd  shell  so  that, 


P  =  cos  1  jl  -  (1  2)1  +  r)/(rMode  +  r)  J  } 


(12) 


By  inspection,  this  gives  an  angle  of  60°  for  each  i-i-i  triplet,  in  agreement  with  angles  of  61.0°, 
60.3°  and  59.0°  from  CRMC  for  Ca-,  Mg-  and  Cu  triplets.  The  Ca-Mg-Ca  (61.1°)  and  Ca-Cu-Ca 
(61.6°)  angles  are  both  slightly  larger  than  60°.  While  this  trend  is  matched  using  Equation  12 
for  these  triplets,  the  estimated  angles  are  somewhat  larger  (63.1°  and  67.6°,  respectively). 
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A  third  peak  in  the  second  coordination  shell  may  be  expected  near  90°  by  summing  the  1st  and 
2nd  peaks.  While  a  local  maxima  is  shown  in  all  five  bond  angle  distributions,  it  is  a  very  shallow 
and  broad  distribution  of  angles.  The  many  configurations  for  intervening  atoms  make  it  difficult 
to  perform  a  simple  topological  analysis  of  the  included  angles. 

3.5  Medium  Range  Order 

The  neutron  diffraction  pre-peak  at  Q  ~  1.2  A'1  of  Ca-Mg-Cu  glasses  (Figure  1)  corresponds  to 
certain  medium-range  correlations  in  real  space  [22].  Generally,  MRO  is  seen  in  solute-lean 
glasses  [15,59,60,61,62]  and  is  explained  in  terms  of  a  sublattice-like  pattern  formed  by  solutes 
[42]  and  significantly  higher  scattering  amplitude  for  solute  than  for  solvent  atoms  [22]. 
However,  similar  to  this  work,  a  pre-peak  has  also  been  seen  in  concentrated  glasses,  such  as 
CessAUs  [63]  and  CaeoMgxZn^-x  [23],  Our  analysis,  given  below  for  the  Ca45Mg25Cu3o  metallic 
glass,  shows  that  the  MRO  in  Ca-Mg-Cu  metal  glasses  is  produced  by  a  face  center  cubic  (FCC) 
-  like  local  ordering  of  Cu-centered  clusters,  and  the  pre-peak  results  from  this  FCC-like 
ordering  combined  with  the  strong  neutron  scattering  from  Cu. 

Figure  12  shows  six  CRMC-simulated  partial  structure  factors,  Aij{Q),  for  the  Ca45Mg25Cu3o 
metallic  glass.  The  dashed  vertical  line  in  these  figures  corresponds  to  the  position  of  the  pre¬ 
peak  maximum  ( Q  =  1.25  A'1)  on  the  respective  total  S(Q )  (see  Figure  la).  It  can  be  clearly  seen 
that  this  pre-peak  is  originated  from  strong  Cu-Cu  and  Mg-Cu  correlations  (Figure  12e  and 
Figure  12f),  as  well  as  a  weak  Mg-Mg  correlation  (Figure  12d).  No  peak  intensity  is  seen  in  this 
Q  range  for  Ca-Ca,  Ca-Mg  and  Ca-Cu  correlations.  This  result  indicates  that,  similar  to  solute- 
lean  metallic  glasses,  MRO  in  Ca-Mg-Cu  metal  glasses  is  caused  by  solute-solute  (mainly  Cu-Cu 
and  Mg-Cu)  interactions. 

Following  a  recent  model  [43],  we  assume  that  MRO  is  caused  by  local  arrangement  of  Cu- 
centered  clusters  in  an  FCC  motif,  as  shown  schematically  in  Figure  13.  This  is  supported  not 
only  by  diffraction  results,  but  also  by  the  triplet  angle  distributions  for  the  2nd  coordination 
shell,  which  show  a  surprising  preference  for  angles  of  60°  and  90°  irrespective  of  composition. 
An  FCC  cell  consists  of  4  Cu-centered  clusters  with  4  octahedral  (P)  and  8  tetrahedral  (y) 
interstitial  sites,  which  tentatively  can  also  be  occupied  by  alloying  elements.  On  average,  each 
Cu-centered  cluster  consists  of  1 1 .4  atoms,  so  that  the  number  of  atoms  per  unit  cell  is  estimated 
to  be  between  45.6  (no  occupied  interstitial  sites)  and  57.6  (all  12  sites  are  occupied).  To  satisfy 
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the  alloy  density  ( pQ  =  0.0373  A'3)  the  cell  volume,  Vc,  should  therefore  be  between  1223  A3  and 
1544  A3,  or  the  cluster  unit  cell  parameter,  ac,  should  be  between  10.7  A  and  1 1.6  A.  This  results 
in  the  closest  average  distance  between  the  Cu-centered  clusters,  rCu-cuCluster  (=«c  N 2) to  be 
between  7.6  A  and  8.2  A.  This  distance  corresponds  to  the  third  peak  on  the  Cu-Cu  PRDF, 
located  at  r  =  7.8  A  (see  Figure  6f).  Assigning  this  value  to  rcu-cuCluster,  the  parameter  ac  is  refined 
as  ac  =  1 1.0  A,  which,  at  the  given  alloy  density,  gives  49.6  atoms  per  super-cell.  This  means  that 
only  4  of  the  12  interstitial  sites  (i.e.  -33%)  per  unit  cell  are  occupied.  Conducting  similar 
analysis  for  the  Mg-centered  clusters  gives  the  average  distance  between  these  clusters,  rMg- 
MgCiustei,  in  the  range  of  8. 1  A  to  8.6  A.  No  peak  is  present  on  the  Mg-Mg  PRDF  in  this  r  range 
(see  Figure  6d).  We  conclude  that  Mg-Mg  clusters  are  not  arranged  in  an  FCC  motif. 

While  the  third  peak  on  the  Cu-Cu  PRDF  can  be  assigned  to  rcu-cuCluster,  which  is  the  distance 
between  the  Cu-centered  clusters  in  the  {0,0,0}  and  { 1/2, 1/2,0}  positions  of  the  FCC  super-cell, 
the  wide  second  peak  on  the  Cu-Cu  PRDF  can  be  assigned  to  the  distances  between  the  Cu  in  the 
center  of  the  cluster  and  Cu  located  in  the  nearest  tetrahedral  and/or  octahedral  sites,  as  well  as  to 
the  distance  between  nearest  tetrahedral  and  octahedral  sites.  The  small  fourth  Cu-Cu  PRDF 
peak  located  at  r  =  9.2  A  can  be  assigned  to  the  distance  between  a  tetrahedral  site  and  second 
nearest  octahedral  site  or  second  nearest  Cu-centered  cluster.  From  crystallography,  the  center  of 
a  tetrahedral  site  in  an  FCC  lattice  is  located  at  a  distance  rTi  =  0.433ac  of  the  nearest  cell  atom 
(Cu-centered  cluster  in  our  case)  and  rT2  =  0.83ac  of  the  second  nearest  cell  atom,  while  the 
center  of  the  octahedral  site  is  located  at  the  distance  ro  =  0.5nc  of  the  nearest  cell  atom.  The  two 
characteristic  distances  between  the  octahedral  and  tetrahedral  sites  are  ron  ~  0.433oc  and  ron  = 
0.83ac.  Therefore,  at  ac=  1 1.0  A,  rTi  =  ron  =  4.8  A  and  ro  =  5.5  A  correspond  to  the  second 
PRDF  peak,  and  rT2  =  r0n  =  9. 13  A  likely  corresponds  to  the  fourth  Cu-Cu  PRDF  peak.  Figure 
14  shows  good  correlation  between  the  Cu-Cu  distances  in  the  FCC  super-lattice,  formed  of  the 
Cu-centered  clusters,  and  the  peak  maxima  on  the  Cu-Cu  partial  radial  distribution  function.  It 
can  also  be  seen  from  this  figure  that  the  Cu-Cu  MRO  correlations  are  only  present  within  one 
FCC  super-cell  and  they  disappear  at  distances  higher  than  ac. 

The  interstitial  sites  can  tentatively  be  occupied  not  only  by  Cu  but  also  Mg  and/or  Ca.  Indeed,  a 
detailed  analysis,  similar  to  the  presented  above,  has  led  to  a  conclusion  that  the  second  peak  in 
the  Mg-Cu  PRDF  can  be  explained  by  the  presence  of  Mg  inside  the  tetrahedral  and  octahedral 
sites.  However,  the  second  peak  in  the  Ca-Cu  PRDF  can  be  explained  by  the  presence  of  Ca 
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inside  octahedral  sites  only.  There  are  no  peaks  in  the  Ca-Cu  PRDF  corresponding  to  the 
tetrahedral  positions.  This  observation  can  be  easily  explained  by  a  large  size  of  Ca  atoms,  which 
would  not  fit  the  small  tetrahedral  sites  [64]. 

3.6  Topological  Description  of  Ca-Mg-Cu  Atomic  Structures 

The  present  results  give  the  most  detailed  experimental  picture  of  a  ternary  metallic  glass  to  date, 
and  so  it  is  worth  analyzing  the  current  findings  with  the  efficient  cluster  packing  (ECP)  model 
for  the  atomic  structure  of  metallic  glasses  [42,43,65].  In  binary  glasses,  preferred  radius  ratios 
and  strong  chemical  interactions  between  unlike  atoms  produce  efficiently  packed  solute- 
centered  clusters  that  are  arranged  in  space  with  an  FCC-like  structure  over  a  length  scale  of  a 
few  cluster  diameters.  In  fonning  these  clusters,  solutes  first  tie  up  the  other  atoms  in  the 
structure,  and  higher  solute  concentrations  give  structures  where  the  cluster-octahedral  and 
cluster-tetrahedral  interstices  are  progressively  filled.  At  even  higher  solute  concentrations, 
solutes  replace  solvent  atoms  in  the  cluster  1st  shell.  Structures  where  all  of  the  cluster-interstitial 
sites  and  about  1/3  of  the  1st  shell  are  occupied  by  solvent  atoms  generally  give  the  most  stable 
binary  glasses.  A  range  of  experimental  and  computational  data  supports  this  model. 

It  was  earlier  suggested  that  the  largest  solute  in  a  multi-component  glass  would  be  the  primary 
structure-forming  solute,  since  larger  solutes  are  topologically  more  potent  [65].  However,  the 
present  work  supports  and  earlier  suggestion  [19,  66]  that  elements  with  the  strongest  chemical 
interaction  with  solvent  atoms  will  fonn  the  structure-forming  clusters. 

While  much  of  the  earlier  analyses  depend  upon  the  nominal  radius  ratio  between  atoms,  it  has 
more  recently  been  shown  that  the  weighted  average  radius  of  atoms  in  the  first  coordination 
shell  gives  a  good  indication  of  GFA  [67].  This  effective  radius  ratio,  f ,  correlates  strongly  with 
good  GFA  and  gives  an  equivalent  structural  description  regarding  local  partial  coordination 
numbers  [64].  This  concept  is  used  here  to  successfully  model  triplet  angles.  Although  the 
nominal  radius  ratio  for  Cu  and  Ca  is  about  0.63  supporting  a  CN=8,  the  effective  radius  ratio  is 
about  0.72.  This  suggests  efficiently  packed  clusters  with  CN=9,  which  is  less  than  the  value  of 
~10  in  the  present  work.  The  cut-off  distance  in  the  present  work  is  about  30%  longer  than  rmocje 
for  Ca  and  Cu  and  50%  larger  for  Mg,  which  may  contribute  to  this  discrepancy. 

The  present  data  show  different  site  filling  rules  for  ternary  glasses.  Inter-cluster  sites  are 
essentially  vacant  in  Ca-Mg-Cu  glasses  with  Cu  concentrations  of  15%  and  25%  (Table  14),  but 
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Cu-Cu  contact  occurs  in  all  glasses  (Table  2),  suggesting  that  some  Cu  atoms  occupy  sites  in  the 
1st  coordination  shell  before  inter-cluster  sites  are  filled.  Topological  analysis  suggests  that 
sufficient  Cu  exists  to  begin  filling  some  cluster-interstitial  sites  at  15%  Cu,  and  the  absence  of 
this  structural  feature  is  consistent  with  the  presence  of  Cu  atoms  in  the  1st  shell  of  Cu-centered 
clusters. 

Solute-rich  binary  metallic  glasses  with  solute  atoms  in  the  1st  coordination  shell  of  solute- 
centered  clusters  are  expected  to  have  4  cluster-interstitial  atoms  per  structure-forming  cluster 
[42,43,65],  or  16  cluster-interstitial  atoms  per  FCC  unit  cell.  However,  recent  phenomenological 
analyses  of  multi-component  BMGs  show  that  they  contain  only  1  or  3  solute  atoms,  or  ‘glue 
atoms’,  per  cluster  [67,68].  In  agreement  with  this,  the  present  work  shows  1  cluster-interstitial 
site  per  Cu-centered  cluster  in  Ca45Mg25Cu3o  and  1.5  sites  per  cluster  for  Ca40Mg25Cu35. 

With  this  as  background,  the  present  data  gives  the  following  topological  picture  of  Ca-Mg-Cu 
structures.  Cu-centered  clusters  in  Ca60Mg25Cui5  have  a  total  CN  of  9.8,  so  that  the  bulk 
concentration  of  Cu  is  achieved  with  1 .7  Cu  atoms  per  cluster.  There  are  no  Cu  atoms  in  cluster- 
interstitial  sites  (Table  14),  so  there  is  1  Cu  atom  at  the  center  of  the  cluster  and  0.7  Cu  atoms  per 
cluster  in  the  1st  coordination  shell.  This  agrees  very  well  with  CNcu-Cu  of  0.6  (Table  2).  The 
composition  of  Cu-centered  clusters  in  Ca6oMg25Cui5  is  nearly  equal  to  the  bulk  composition, 
further  supporting  this  interpretation.  The  total  CN  around  Cu  clusters  is  10.4  in  Ca50Mg25Cu25, 
so  there  are  3.5  Cu  atoms  per  structure-forming  cluster.  Data  in  Table  14  suggest  there  are  no 
atoms  in  cluster-interstitial  sites,  so  that  there  would  be  1  Cu  atom  in  the  cluster  center  and  2.5 
Cu  atoms  in  the  1st  coordination  shell.  However,  this  gives  a  partial  CN  that  is  larger  than  the 
value  in  Table  2,  and  the  composition  of  Cu-centered  clusters  is  noticeably  deficient  in  Cu.  Since 
the  data  in  Table  14  has  the  largest  uncertainty,  this  structural  interpretation  is  unlikely.  The  bulk 
composition  is  achieved  when  0.7  Cu  atoms  in  cluster-interstitial  sites  are  added  to  the  Cu- 
centered  cluster  composition,  so  that  an  alternate  interpretation  is  that  there  is  1  Cu  atom  at  the 
center  of  the  clusters,  0.7  Cu  atoms  in  cluster-interstitial  sites,  and  1.8  Cu  atoms  in  the  1st 
coordination  shell  of  the  Cu-centered  clusters.  This  structure  gives  a  partial  CNcu-Cu  that  is 
consistent  with  the  data  in  Table  2  and  with  the  observation  of  a  peak  in  the  Cu-Cu  g(r)  near  7.8 
A.  The  total  Cu  CN  is  10.4  in  Ca45Mg25Cu3o,  so  there  are  4.5  Cu  atoms  per  cluster.  The 
composition  of  Cu-centered  clusters  is  deficient  in  Cu,  so  that  an  additional  0.7  Cu  atoms  per  are 
needed  to  match  the  bulk  composition.  This  suggests  1  Cu  atom  at  the  cluster  center,  0.7  Cu 
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atoms  at  cluster-interstitial  sites  (Table  14),  and  2.8  Cu  atoms  in  the  1st  coordination  shell.  This  is 
consistent  with  the  data  in  Table  14,  with  bulk  composition  and  partial  CNs,  and  with  the  Cu-Cu 
peak  in  g(r)  near  7.8  A.  Finally,  there  are  5.5  Cu  atoms  per  structure-forming  cluster  in 
Ca4oMg25Cu35.  To  match  bulk  composition,  1.8  Cu  atoms  must  occupy  cluster-interstitial  sites. 
Thus,  there  is  1  Cu  atom  at  the  cluster  centers,  1.8  Cu  atoms  in  cluster-interstitial  sites,  and  2.7 
Cu  atoms  in  the  1st  coordination  shell.  This  structure  is  consistent  with  the  number  of  occupied 
interstitial  sites  from  Table  14  (1.6±0.3),  with  bulk  composition  and  Cu-Cu  partial  CNs  (Table 
2),  and  with  the  Cu-Cu  g(r). 

While  bulk  compositions  can  be  matched  by  adding  Cu  atoms  to  cluster-interstitial  sites,  the  Mg 
concentrations  are  still  slightly  deficient.  The  present  work  suggests  that  some  Mg  is  in  cluster- 
interstitial  sites,  consistent  with  this  observation. 

CONCLUSIONS 

1.  The  atomic  structures  of  four  Ca4o+xMg25Cu35_x  (X  =  0,  5,  10,  and  20  at.%)  ternary  metallic 
glasses  have  been  determined  using  a  synergistic  combination  of  neutron  diffraction,  quantum 
molecular  dymamics  (QMD)  and  constrained  reverse  Monte  Carlo  (CRMC)  modeling.  All  six 
partial  radial  distribution  functions  (PRDFs),  gij(r),  have  been  identified  for  each  alloy.  The 
S(Q)  curves  damp  very  quickly  with  Q  and  no  oscillations  are  apparent  beyond  ~15  A'1, 
indicating  that  the  glasses  have  a  wide  range  of  interatomic  distances.  An  increase  in  Cu 
content  shortens  the  average  interatomic  distance  and  narrows  the  distribution  of  interatomic 
distances,  shifting  the  first  sharp  diffraction  peak  to  higher  Q  and  broadening  the  peak.  A  pre¬ 
peak  at  Q  ~  1.2  A'1  in  the  S(Q)  curves  indicates  the  presence  of  medium  range  order  (MRO)  in 
these  amorphous  alloys. 

2.  The  first  PRDF  peaks,  which  correspond  to  interatomic  bond  distances  in  the  1st  coordination 
shell,  are  asymmetric,  and  the  mode  bond  distances  are  always  smaller  than  the  respective 
mean  bond  distances.  The  difference  is  the  smallest  for  Ca-Ca  pairs  and  increases  in  the  order 
of  Ca-Mg,  Mg-Mg,  Ca-Cu,  Cu-Cu  and  Mg-Cu.  This  indicates  that  Ca-Ca  atom  pairs  are  most 
often  compressed  but  often  show  ‘hard-sphere’  contact,  and  that  Mg-Mg  and  Cu-Cu 
displacements  most  commonly  show  ‘hard-sphere’  contact  but  often  have  extended  bond 
distances  within  the  first  shell.  Atomic  separations  for  unlike  atom  pairs  are  slightly  shorter 
than  the  sum  of  atomic  radii,  suggesting  that  the  atoms  are  slightly  non-spherical. 
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3.  The  nearest-neighbor  bond  lengths  are  shorter  than  those  in  competing  crystals.  A  noticeable 
shortening  of  Ca-Cu  and  Mg-Cu  bond  distances  indicate  strong  interactions  between  these 
atom  pairs.  It  is  suggested  that  the  bond  shortening  is  enabled  by  the  absence  of  long-range 
atomic  order,  which  lowers  the  free  energy  of  metallic  glasses  and  increases  GFA. 

4.  Pronounced  chemical  short  range  ordering  (CSRO)  near  Cu  atoms,  chemical  short  range 
clustering  (CSRC)  near  Ca  atoms  and  a  neutral  environment  near  Mg  atoms  are  shown. 
Increasing  the  Cu  concentration  from  15  to  35  atomic  percent  increases  the  total  coordination 
number  around  Ca,  CNca,  from  13.6  to  15.0,  while  CNMg  and  CNcu  remain  unchanged  at  those 
-12.5  and  -10.4,  respectively.  The  partial  coordination  numbers  depend  nearly  linearly  on 
alloy  composition,  so  that  the  number  of  Ca  atoms  decrease,  the  number  of  Cu  atoms  increase 
and  the  number  of  Mg  atoms  are  almost  constant  with  increasing  Cu  concentration. 

5.  Voronoi  tessellation  and  three-body  correlations  show  that  many  types  of  coordination 
polyhedra  are  present,  but  the  most  common  are  (0,2, 8, 4)  Ca-centered,  (0,2, 8, 2)  Mg-centered 
and  (0,3, 6,0)  and  (0,2,8, 1)  Cu-centered.  The  fractions  and  distributions  of  these  clusters  depend 
on  alloy  composition.  Polytetrahedral-type  clusters  and  five-coordinated  vertices  dominate  in 
the  amorphous  structures,  which  indicate  that  tetrahedra  and  pentagonal  bi-pyramids  are  the 
main  building  blocks  in  these  amorphous  alloys.  The  GFA  of  the  Ca-Mg-Cu  alloys  increases 
with  a  decreasing  fraction  of  Cu-centered  polytetrahedral  clusters.  The  local  atom  packing 
fraction  near  a  given  atom  increases  with  atomic  radius  and  with  an  increase  in  Cu 
concentration.  Global  atom  packing  fractions  of  0.68  to  0.73  are  much  larger  than  for  dense 
random  packed  structures,  and  are  comparable  to  packing  in  crystalline  structures. 

6.  Medium  range  order  is  present  in  Ca-Mg-Cu  amorphous  alloys  and  is  the  result  of  strong 
tendency  of  Cu  for  CSRO.  Detailed  analysis  shows  that  the  amorphous  structure  of  these  alloys 
is  described  by  close-packing  of  Cu-centered  clusters  that  follows  a  face-centered  cubic  motif 
over  a  length  scale  of -10  A,  or  about  3  cluster  diameters.  This  gives  a  characteristic  medium 
range  order  pre-peak  at  Q  -  1.2  A  in  the  total  scattering  structure  factors  of  these  alloys. 

7.  A  topological  description  of  Ca-Mg-Cu  atomic  structures  shows  that  the  effective  radius  of 
atoms  in  the  1st  coordination  shell  of  Cu  atoms  enables  efficient  atomic  packing.  These  clusters 
match  the  bulk  composition  for  Ca6oMg2sCui5,  but  are  Cu-lean  for  all  other  compositions.  The 
remaining  structures  have  between  0.7  and  1.8  Cu  atoms  in  cluster-interstitial  sites.  This  allows 
the  bulk  composition  to  be  matched  and  is  necessary  for  agreement  with  the  Cu-Cu  g(r)  plots. 
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Cu  atoms  occupy  between  0.7  (for  Ca6oMg25Cui5)  and  2.8  (for  Ca45Mg25Cu3o)  sites  in  the  first 
coordination  shells  of  Cu-centered  clusters. 
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TABLES 


Table  1.  Density  (in  g/cnr  and  atoms/A3)  of  Ca-Mg-Cu  amorphous  alloys  studied  in  this  work. 


Alloy 

Density 

g/cnr 

Atoms/A3 

Ca4oMg25Cu35 

2.936  ±0.007 

0.03987 

Ca45Mg25Cu3o 

2.673  ±0.005 

0.03728 

Ca5oMg25Cu25 

2.439  ±0.003 

0.03497 

Ca60Mg25Cui5 

2.039  ±0.003 

0.03097 

Table  2.  Partial  and  total  coordination  numbers  around  Ca,  Mg  and  Cu  atoms  in  the  Ca-Mg-Cu 
amorphous  alloys,  in  accord  to  QMD  and  RMC  simulations  and  Voronoi  (Vor)  tessellation. 


Ca. 

wMg25CU35 

Ca45Mg25Cu3o 

Ca. 

«Mg2sCU25 

Ca, 

SoMg25CUi5 

^iN 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

QMD 

CRMC 

Vor 

Ca-Ca 

7.0 

6.9 

7.0 

7.8 

7.1 

7.1 

8.6 

7.9 

7.9 

9.2 

8.6 

8.6 

Ca-Mg 

4.3 

3.7 

3.6 

4.1 

3.7 

3.8 

3.8 

3.4 

3.4 

3.5 

3.3 

3.4 

Ca-Cu 

4.9 

4.4 

4.4 

3.9 

3.8 

3.8 

3.2 

3.1 

3.1 

1.6 

1.7 

1.7 

Total 

16.2 

15.0 

15.0 

15.8 

14.6 

14.7 

15.6 

14.4 

14.4 

14.3 

13.6 

13.7 

Mg-Ca 

6.9 

5.9 

5.8 

7.3 

6.6 

6.9 

7.7 

6.8 

6.9 

8.5 

8.0 

8.0 

Mg-Mg 

2.6 

3.1 

3.0 

2.3 

2.8 

2.8 

2.1 

2.7 

2.8 

1.8 

2.7 

2.6 

Mg-Cu 

3.2 

3.7 

4.0 

2.5 

3.4 

3.3 

2.1 

3.0 

2.8 

1.1 

1.3 

1.6 

Total 

12.8 

12.7 

12.8 

12.2 

12.8 

13.0 

11.8 

12.5 

12.5 

11.3 

12.0 

12.2 

Cu-Ca 

5.6 

5.0 

5.1 

5.9 

5.6 

5.7 

6.4 

6.2 

6.2 

6.4 

6.8 

6.7 

Cu-Mg 

2.3 

2.7 

2.9 

2.1 

2.8 

2.7 

2.1 

3.0 

3.0 

1.8 

2.4 

2.6 

Cu-Cu 

1.6 

2.6 

2.1 

1.3 

1.9 

2.0 

1.0 

1.2 

1.5 

0.6 

0.6 

0.8 

Total 

9.4 

10.3 

10.1 

9.3 

10.4 

10.4 

9.5 

10.4 

10.7 

8.8 

9.8 

10.1 

24 

Approved  for  public  release;  distribution  unlimited 


Table  3.  The  most  frequent  (mode)  distances  (in  Angstroms)  between  pair  atoms  in  Ca-Mg-Cu 
amorphous  alloys,  in  accord  to  QMD  and  CRMC  simulations. 


Alloy 

Method 

Ca-Ca 

Ca-Mg 

Ca-Cu 

Mg-Mg 

Mg-Cu 

Cu-Cu 

Ca4oMg25CU35 

QMD 

3.75 

3.42 

3.02 

3.14 

2.71 

2.47 

CRMC 

3.76 

3.45 

3.01 

3.14 

2.61 

2.49 

Ca45Mg25Cu30 

QMD 

3.75 

3.41 

3.03 

3.13 

2.71 

2.48 

CRMC 

3.74 

3.43 

3.04 

3.12 

2.62 

2.49 

Ca5oMg25CU25 

QMD 

3.75 

3.41 

3.07 

3.13 

2.75 

2.48 

CRMC 

3.75 

3.43 

3.04 

3.13 

2.73 

2.46 

Ca60Mg25Cui5 

QMD 

3.78 

3.46 

3.03 

3.08 

2.72 

2.51 

CRMC 

3.78 

3.47 

3.08 

3.12 

2.70 

2.50 

Average 

CRMC 

3.76 

3.45 

3.04 

3.12 

2.67 

2.49 

Table  4.  Weighted  average  (mean)  distances  (in  Angstroms)  between  pair  atoms  in  the  first  shell 
in  Ca-Mg-Cu  amorphous  alloys,  in  accord  to  QMD  and  CRMC  simulations. 

Alloy _ Ca-Ca  Ca-Mg  Ca-Cu  Mg-Mg  Mg-Cu  Cu-Cu 

„  .  .  „  QMD 

Ca40Mg25Cu35  (CRJVIC 

„  ,  ,  „  QMD 

Ca45Mg25CU30 


Ca5oMg25CU25 

QMD 

CRMC 

4.01 

3.89 

3.62 

3.59 

3.25 

3.33 

3.24 

3.30 

2.89 

3.15 

2.58 

2.63 

Ca60Mg25Cui5 

QMD 

4.00 

3.64 

3.19 

3.29 

2.87 

2.60 

CRMC 

3.89 

3.66 

3.38 

3.46 

2.75 

2.60 

Average 

CRMC 

3.88 

3.60 

3.29 

3.37 

3.03 

2.71 

Table  5.  The  first  shell  cut-off  distances  between  pair  atoms  (in  Angstroms)  in  Ca-Mg-Cu 

amorphous  alloys,  in  accord  to  QMD  and  CRMC  simulations. 

Alloy 

Ca-Ca 

Ca-Mg 

Ca-Cu  1 

Mg-Mg 

Mg-Cu 

Cu-Cu 

Ca4oMg25Cu35 

QMD 

4.88 

4.63 

4.26 

4.22 

3.58 

3.25 

CRMC 

4.74 

4.45 

4.11 

4.29 

4.04 

3.75 

Ca45Mg25Cu30 

QMD 

4.91 

4.59 

4.26 

4.09 

3.58 

3.34 

CRMC 

4.73 

4.63 

4.31 

4.26 

4.34 

3.54 

Ca5oMg25Cu25 

QMD 

4.99 

4.58 

4.40 

3.97 

3.84 

3.42 

CRMC 

4.82 

4.54 

4.44 

4.20 

4.30 

3.32 

Ca60Mg25Cui5 

QMD 

4.94 

4.58 

4.13 

4.05 

3.89 

3.31 

CRMC 

4.84 

4.70 

4.41 

4.46 

3.92 

3.22 

3.95 

3.64 

3.21 

3.31 

2.83 

2.59 

3.89 

3.55 

3.19 

3.34 

3.03 

2.85 

3.98 

3.63 

3.22 

3.28 

2.84 

2.59 

3.85 

3.61 

3.27 

3.37 

3.19 

2.75 
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Table  6.  Voronoi  cell  volume,  void  volume  and  packing  fraction  of  Ca,  Mg,  and  Cu  atoms  in  the 
amorphous  structure  of  Ca-Mg-Cu  alloys.  The  packing  fractions  of  the  alloys  are  also  given  in 
the  last  column. 


Table  9.  Metallic,  rm,  [69]  and  covalent,  rc,  [70]  crystalline  bond  distances  between  Ca,  Mg  and 
Cu  atom  pairs.  These  distances  are  compared  to  rMode  and  rMean  from  the  CRMC  PRDFs.  rMode 
and  rMean  are  the  average  values  for  the  four  studied  alloys. 


Ca-Ca 

Ca-Mg 

Ca-Cu 

Mg-Mg 

Mg-Cu 

Cu-Cu 

Tm  (A) 

3.94 

3.57 

3.25 

3.20 

2.88 

2.56 

rc  (A) 

3.52 

3.17 

3.08 

2.82 

2.73 

2.64 

I*Mode  C^-) 

3.76 

3.45 

3.04 

3.13 

2.67 

2.49 

I*Mean  (-A.) 

3.88 

3.60 

3.29 

3.37 

3.03 

2.71 

100  /o(rMean/l*Mode"  1 ) 

3.2 

4.3 

8.2 

7.7 

13.5 

OO 

OO 

1 00%(rM/rMode- 1 ) 

4.8 

3.5 

6.9 

2.2 

7.9 

2.8 

100%(rc/rMode-l) 

-6.4 

-8.1 

1.3 

-9.9 

2.2 

6.0 

100%(rM/rMean-l) 

1.5 

-0.8 

-1.2 

-5.0 

-5.0 
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Table  10.  Linear  fits  to  the  dependence  on  Cu  content  (at.%)  of  coordination  number  data 
(Figure  1 1).  Rf  is  the  coefficient  of  determination  of  the  respective  linear  regression. 


Atomic  pair 

Gradient 

Intercept 

R2 

Ca-Ca 

-0.090 

9.98 

0.963 

Ca-Mg 

0.022 

2.95 

0.884 

Ca-Cu 

0.125 

0.00 

0.992 

Mg-Ca 

-0.101 

9.48 

0.978 

Mg-Mg 

0.018 

2.36 

0.639 

Mg-Cu 

0.110 

0.00 

0.956 

Cu-Ca 

-0.089 

8.24 

0.966 

Cu-Mg 

0.0154 

2.32 

0.278 

Cu-Cu 

0.099 

-1.02 

0.951 

Table  11.  Short-range  order  parameters,  a,,,  for  Ca-,  Mg-  and  Cu-  centered  clusters  in  the  Ca- 
Mg-Cu  amorphous  alloys. 


aij 

Ca-Ca 

Ca-Mg 

Ca-Cu 

Mg-Ca 

Mg-Mg 

Mg-Cu 

Cu-Ca 

Cu-Mg 

Cu-Cu 

Ca4oMg25Cu35 

-0.16 

0.03 

0.16 

-0.15 

0.00 

0.17 

-0.21 

-0.05 

0.28 

Ca45Mg25Cu3o 

-0.07 

-0.04 

— 

-0.17 

0.15 

0.14 

-0.21 

-0.08 

0.38 

Ca5oMg25Cu25 

-0.10 

0.05 

mam 

-0.09 

0.14 

0.05 

-0.15 

-0.12 

0.42 

Ca6oMg25Cui5 

-0.05 

0.02 

0.17 

-0.11 

0.11 

0.25 

-0.17 

0.07 

0.56 

Table  12.  The  short-range  order  parameters,  a!(Jk),  for  Ca-,  Mg-  and  Cu-  centered  clusters  in  the 


Ca-Mg-Cu  amorphous  alloys. 


Ca 

Mg 

Cu 

Ca40Mg25Cu35 

0.11 

0.00 

-0.15 

Ca45Mg25Cu3o 

0.06 

-0.05 

-0.16 

Ca5oMg25Cu25 

0.10 

-0.05 

-0.14 

Ca60Mg25Cui5 

0.08 

-0.04 

-0.10 
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Table  13.  Fractions  of  Ca-,  Mg-  and  Cu-centered  Kasper  polyhedra,  distorted  Kasper  polyhedra 
and  sum  of  both. 


Distorted  Kasper 

Kasper  Polyhedra 

Polyhedra 

Sum 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Ca 

Mg 

Cu 

Ca40Mg25CU35 

0.059 

0.142 

0.323 

0.408 

0.467 

0.275 

0.467 

0.610 

0.598 

Ca45Mg25Cu3o 

0.079 

0.152 

0.323 

0.424 

0.402 

0.269 

0.503 

0.555 

0.592 

Ca5oMg25CU25 

0.102 

0.192 

0.308 

0.454 

0.425 

0.260 

0.556 

0.617 

0.568 

Ca60Mg25Cui5 

0.121 

0.190 

0.338 

0.498 

0.407 

0.283 

0.619 

0.597 

0.621 

Table  14.  The  lattice  parameter,  ac,  the  number  of  atoms  per  unit  cell,  Nc,  and  the  number  of 
occupied  interstitial  sites  per  unit  cell,  /V,-  in  the  face  centered  cubic  lattice  formed  of  the  Cu- 


centered  clusters  in  Ca-Mg-Cu  amorphous  alloys. 


Alloy 

Ca4oMg25Cu35 

Ca45Mg25Cu3o 

Ca50Mg25CU25 

Ca4()Mg2sCu  is 

ac  (A) 

10.9  ±0.1 

1 1.0  ±  0.1 

10.9  ±0.1 

11.2  ±0.2 

Nc 

51.6  ±  1.5 

49.6  ±  1.5 

45.3  ±1.3 

43.6  ±2.3 

N, 

6.4  ±  1.5 

4.0  ±  1.5 

0.0  ±  1.3 

1.2  ±2.3 
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FIGURE  CAPTIONS 


Figure  1 .  (a)  Experimental  total  neutron  scattering  structure  factors  for  Ca6o-xMg25Cui5+x 
metallic  glasses.  The  curves  are  separated  by  shifting  along  the  vertical  axis  for  better  view,  (b) 
Dependence  on  the  concentration  of  copper  of  the  positions  of  the  start  (1st  Minimum), 
maximum  (Peak)  and  end  (2nd  minimum)  of  the  first  peak,  as  well  as  the  peak  width. 

Figure  2.  (a)  Total  radial  distribution  functions  for  Ca6o-xMg2sCui5+x  metallic  glasses.  The 
curves  are  separated  by  shifting  along  the  vertical  axis  for  better  view,  (b)  Dependence  on  the 
concentration  of  copper  of  the  positions  of  the  first,  second  and  third  peaks  and  the  width  of  the 
first  and  second  peaks. 

Figure  3.  QMD-simulated  (a)  Cu-Cu,  (b)  Mg-Cu,  (c)Mg-Mg,  (d)  Ca-Cu,  (e)  Ca-Mg  and  (f)  Ca- 
Ca  PRDFs  for  Ca4oMg25Cu35,  Ca50Mg25Cu25  and  Ca6oMg25Cui5  metallic  glasses. 

Figure  4.  Experimental  (dark  solid  lines)  and  QMD-simulated  (red  dashed  lines)  RDFs  for  Ca6o- 
xMg25 Cu 1 5+x  metallic  glasses.  The  curves  are  separated  by  0.3-point  shifting  along  the  vertical 
axis  for  better  view. 

Figure  5.  Experimental  (dark  solid  lines)  and  CRMC-simulated  (red  dashed  lines)  (a)  total 
structure  factors  and  (b)  RDFs  for  Ca6o-xMg25Cui5+x  metallic  glasses.  The  curves  are  separated 
by  shifting  along  the  vertical  axis  for  better  view. 

Figure  6.  CRMC-simulated  (a)  Ca-Ca,  (b)  Ca-Mg,  (c)  Ca-Cu,  (d)  Mg-Mg,  (e)  Mg-Cu  and  (f)  Cu- 
Cu  partial  radial  distribution  functions,  gy(r)  for  Ca6o-xMg25Cui5+x  metallic  glasses. 

Figure  7.  Distribution  of  different  types  of  (a)  Ca-centered  (b)  Mg-centered  and  (c)  Cu-centered 
clusters  in  Ca6o-xMg25Cui5+x  amorphous  alloys. 

Figure  8.  Fractions  of  3-,  4-,  5-  and  6-  coordinated  vertices  in  Ca-,  Mg-  and  Cu-  centered  clusters 
in  (a)  Ca4oMg25Cu35,  (b)  Ca45Mg25Cu3o,  (c)  Ca5oMg25Cu25,  and  (d)  Ca6oMg25Cui5  amorphous 
alloys. 

Figure  9.  Ca-Ca-Ca,  Ca-Mg-Ca,  Ca-Cu-Ca,  Mg-Mg-Mg  and  Cu-Cu-Cu  bond  angle  distributions 
in  (a)  the  first  coordination  shell  and  (b)  second  coordination  shell  in  the  Ca45Mg25Cu3o 
amorphous  alloy. 


30 

Approved  for  public  release;  distribution  unlimited 


Figure  10.  Comparison  of  the  mode  atomic  pair  distances  with  (a)  the  respective  mean  distances 
in  the  amorphous  structure  of  the  studied  Ca-Mg-Cu  alloys,  and  with  the  reported  (b)  metallic, 
rM,  and  (c)  covalent,  rc,  bond  distances. 

Figure  1 1 .  Dependence  on  the  Cu  concentration  of  (a)  the  total  coordination  number,  Nm,  where 
M  is  Ca,  Mg,  or  Cu  center  atom,  and  partial  coordination  numbers  (b)  Nm- Ca,  (c)  AVi-Mg,  and  (d) 

Mvi-cu- 

Figure  12.  Partial  structure  factors,  Aij(Q),  for  Ca45Mg25Cu3o  metallic  glass. 

Figure  13.  Schematic  of  a  face  centered  cubic  cell  fonned  by  Cu-centered  clusters  (shown  by 
red).  Positions  of  octahedral  (P)  and  tetrahedral  (y)  interstitial  voids  are  also  shown. 

Figure  14.  Correlation  of  the  Cu-Cu  distances  in  the  FCC  super-lattice  formed  of  the  Cu-centered 
clusters  and  the  intensity  peaks  on  the  Cu-Cu  partial  radial  distribution  function. 
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FIGURES 


Figure  1 .  (a)  Experimental  total  neutron  scattering  structure  factors  for  Ca6o-xMg25Cui5+x 
metallic  glasses.  The  curves  are  separated  by  shifting  along  the  vertical  axis  for  better  view,  (b) 
Dependence  on  the  concentration  of  copper  of  the  positions  of  the  start  (1st  Minimum), 
maximum  (Peak)  and  end  (2nd  minimum)  of  the  first  peak,  as  well  as  the  peak  width. 
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Figure  2.  (a)  Total  radial  distribution  functions  for  Ca6o-xMg2sCui5+x  metallic  glasses.  The 
curves  are  separated  by  shifting  along  the  vertical  axis  for  better  view,  (b)  Dependence  on  the 
concentration  of  copper  of  the  positions  of  the  first,  second  and  third  peaks  and  the  width  of  the 
first  and  second  peaks. 
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Figure  3.  QMD-simulated  (a)  Cu-Cu,  (b)  Mg-Cu,  (c)Mg-Mg,  (d)  Ca-Cu,  (e)  Ca-Mg  and  (f)  Ca- 
Ca  PRDFs  for  Ca4oMg25Cu35,  Ca5oMg25Cu25  and  Ca6oMg2sCui5  metallic  glasses. 
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Figure  4.  Experimental  (dark  solid  lines)  and  QMD-simulated  (red  dashed  lines)  RDFs  for 
Ca6o-xMg25Cui5+x  metallic  glasses.  The  curves  are  separated  by  0.3-point  shifting  along  the 
vertical  axis  for  better  view. 
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Figure  5.  Experimental  (dark  solid  lines)  and  RMC-simulated  (red  dashed  lines)  (a)  total 
structure  factors  and  (b)  RDFs  for  Ca6o-xMg25Cui5+x  metallic  glasses.  The  curves  are  separated 
by  0.3-point  shifting  along  the  vertical  axis  for  better  view. 
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Figure  6.  RMC-simulated  (a)  Ca-Ca,  (b)  Ca-Mg,  (c)  Ca-Cu,  (d)  Mg-Mg,  (e)  Mg-Cu  and  (f)  Cu- 
Cu  partial  radial  distribution  functions,  gy(r)  for  Ca6o-xMg25Cuis+x  metallic  glasses. 
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Figure  7.  Distribution  of  different  types  of  (a)  Ca-centered  (b)  Mg-centered  and  (c)  Cu-centered 
clusters  in  Ca6o-xMg25Cui5+x  amorphous  alloys. 
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Figure  8.  Fractions  of  3-,  4-,  5-  and  6-  coordinated  vertices  in  Ca-,  Mg-  and  Cu-  centered  clusters 
in  (a)  Ca4oMg25Cu35,  (b)  Ca45Mg25Cu3o,  (c)  Ca5oMg25Cu25,  and  (d)  Ca60Mg25Cui5  amorphous 
alloys. 
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Figure  9.  Ca-Ca-Ca,  Ca-Mg-Ca,  Ca-Cu-Ca,  Mg-Mg-Mg  and  Cu-Cu-Cu  bond  angle  distributions 
in  (a)  the  first  coordination  shell  and  (b)  second  coordination  shell  in  the  Ca45Mg2sCu3o 
amorphous  alloy. 
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Figure  10.  Comparison  of  the  mode  atomic  pair  distances  with  (a)  the  respective  mean  distances 
in  the  amorphous  structure  of  the  studied  Ca-Mg-Cu  alloys,  and  with  the  reported  (b)  metallic, 
rM,  and  (c)  covalent,  rc,  bond  distances. 
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Figure  11.  Dependence  on  the  Cu  concentration  of  (a)  the  total  coordination  number,  Nm,  where 
M  is  Ca,  Mg,  or  Cu  center  atom,  and  partial  coordination  numbers  (b)  A'M-ca,  (c)  NM- Mg,  and  (d) 
A'm-Cu- 
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Figure  13.  Schematic  of  a  face  centered  cubic  cell  formed  by  Cu-centered  clusters  (shown  by 
red).  Positions  of  octahedral  (P)  and  tetrahedral  (y)  interstitial  voids  are  also  shown. 
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Figure  14.  Correlation  of  the  Cu-Cu  distances  in  the  FCC  super-lattice,  formed  of  the  Cu- 
centered  clusters,  and  the  intensity  peaks  on  the  Cu-Cu  partial  radial  distribution  function. 
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